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Abstract. We stir vortices into a trapped quasi two-dimensional atomic Bose-Einstein 
condensate by moving three laser stirrers. We apply stirring protocols introduced by Boyland 
\et al. \ ( |2000[ ), that efficiently build in topological chaos in classical fluids and are classified as 
Pseudo-Anosov stirring protocols. These are compared to their inefficient mixing counterparts, 
finite-order stirring protocols. We investigate if inefficient stirring protocols result in a more 
clustered distribution of vortices. The efficiency with which vortices are 'mixed' or distributed in 
a condensate is important for investigating dynamics of continuously forced quantum turbulence 
and the existence of the inverse cascade in turbulent two-dimensional superfluids. 



1. Introduction 

Jupiter's Great Red Spot is a beautiful demonstration of a large spot of vorticity, one of the 
defining features of two dimensional (2D) classical turbulence, arising from the inverse cascade 
process (Sommeria et al. (1988); Marcus (1988)). This process describes the flow of energy 



from small scales, given by the scale at which energy is injected into the system, to larger scales, 
manifesting in the clustering of vortices of like- winding. The concept of an inverse energy cascade 



process was introduced by Onsager (1949), who showed that clusters of like-signed point vortices 



in a large Onsager gas of many point-vortices possess negative temperature. Quantum vortices 
in superfluid He can be realistically modelled as point vortices or vortex lines, due their small 
vortex core size (~ lnm), quantised vorticity and large ratio of intervortex spacing to vortex 
core radius (~ 10 5 — 10 6 ). At large scales greater than the average inter- vortex spacing, three 
dimensional turbulent quantum fluids are also known to exhibit the same classical Kolmogorov 
energy spectrum with scaling k~ 5 ^ 3 (where k is the wavenumber) as classical fluids (see 



Nore 



et al.\ ([1997D; |Maurer fc Tabelingj fll998[); |Stalp et al.\ Ql999[ ) ; |Araki et al\ fl2002P; |Kobayashi & 



Tsubota (2005); Salort et al. (2010); Baggaley k Barenghi (2011) and Baggaley et ol[(2012)). 



For these reasons, one might expect the inverse cascade process to occur in other two dimensional 
quantum fluids, such as Bose-Einstein condensates (BECs), where vortices also have quantised 
circulation. While at large scales turbulent quantum fluids are known to have similarities to 
turbulent classical fluids for three dimensional systems, for two dimensional systems such a 
crossover is yet to be established. There have been some numerical simulations working towards 
this goal in trapped and homogeneous condensates, however, to date there is no conclusive 



demonstration of the inverse cascade process for harmonically trapped atomic BECs (Parker 



& Adams ('. 


2005 


); 


Horng et al. 


( 


2009); 


Numasato et al. 


( 


2010 


White et al. 


( 


2012 


); Bradley & Anderson ( 


2012 


); 


Reeves et al. 



( |2012a|6D ). Unlike in superfluid 



He, atomic BECs of low dimensions can be routinely created and manipulated experimentally 
and so they provide a promising system for such investigations. One important component of 
any two dimensional turbulent set-up is the introduction of vortices into the condensate in a 
controlled manner, such that their distribution is well-mixed. Finding such a well-mixed vortex 
distribution will be the focus of this work. 



Recently White et al. (2012) showed that the trajectory a single laser stirrer takes through the 
condensate determines if the resulting vortices are clustered or more randomly distributed. We 
found that a single laser moving on a circular trajectory creates an initial cluster of like-signed 
vortices. While the extent of clustering was found to increase with stirring, it did not persist 
and decayed after the laser stirrer was turned off. This clustering of vortices manifests due to 
the trajectory of the laser stirrer, and not as a result of an inverse cascade process. Indeed it 
was shown that altering the trajectory of the laser stirrer resulted in less-extensive clustering 
and vortices could even be created in more randomly distributed configurations. In this work we 
employ three laser stirrers and measure the distribution and clustering of vortices by applying 
statistical measures of clustering previously developed in White et al. (2012). In particular, 
we stir the condensate in a manner that is known to efficiently mix classical fluids, known as 
pseudo-Anosov stirring protocols. These are contrasted to finite-order stirring protocols which 
have been shown to be inefficient mixers for classical fluids. 

In this paper we proceed in the following way. Firstly we give an overview of topological 
stirring protocols and then outline the model applied to simulate vortex dynamics in atomic 
BECs. Finally we introduce statistical measures of clustering that are applied to give insight into 
how two quite different stirring protocols influence vortex dynamics in a surprisingly similar way. 
While this is unexpected from the perspective of classical fluids, it is a remarkable demonstration 
of the potential flow of quantum fluids. 



2. Topological stirring protocols 

Topological concepts were introduced to address the question of how to most efficiently mix 
classical fluids. Boyland et al. (2000) determined that for flows with the topology of certain 
braids, given only the topology of the flow, it is possible to deduce the material stretch rate. For 
stokes flow (characterized by Reynolds number, Re <C 1), chaotic advection of fluid particles is 
built into the flow faster by stirring schemes that build in topological chaos. Such protocols are 
classified as pseudo-Anosov stirring protocols and it is established that stirring a fluid in such 
a way results in exponential stretching of material lines. Stirring protocols that do not have 



the correct topology to mix effectively are known as finite-order stirring protocols (Boyland 



et al. 


(2000 


); 


Finn et al. 


(2003 


); 



. In classical fluids, finite-order stirring 
protocols have been shown to stretch material lines linearly and therefore do not mix as efficiently 
as their pseudo-Anosov counterparts. 

The trajectory of the stirring lasers can be represented by braid- diagrams in space-time, with 
braid letters a% denoting the exchange of stirrers labelled % and i + 1 in a clockwise direction. This 
can be depicted in a braid diagram by the strand corresponding to the stirrer labelled i+1 passing 
over the strand showing the ith stirrer. Similarly, an anti-clockwise path of stirrers through the 
condensate is denoted by braid letters &-{ which represents stirrers i + 1 and i exchanging 
position in an anti-clockwise direction. This interchange is shown on a braid-diagram by the 
strand depicting the i + 1 stirrer passing under the ith strand. Figure [I] (A) illustrates the braid 
diagram corresponding to the exchange of stirrers following a protocol that can be classified as 
pseudo-Anosov. 

Each stirring operation, corresponding to a braid letter, can be represented by a braid matrix. 
For the stirring operations corresponding to the a\, G—\ and a 2 braid letters, which we apply 



in this work, the corresponding braid matrices take the form (Boyland et al. (2000); Finn et al. 



(2003)) 
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(1) 

The time-like sequence of braid letters composes a braid word, where the rightmost letter gives 
the first stirring operation and the leftmost letter denotes the most recent stirring operation. 
We concentrate on two stirring protocols, which we label A and B that are represented by the 
braid words 

A = (aia 2 ) n and B = (<7-i<7 2 ) n , (2) 

here n denotes the number of repetitions of braid letter operations. Similarly, the corresponding 
action of the braid matrices takes the form 



A = [SiS 2 ] n and B = [S_iE 2 ] 



(3) 



Stirring protocol A is a finite-order protocol, while protocol B is a pseudo-Anosov stirring 
protocol which builds topological chaos into the flow. The topology of the fluid flow is completely 
characterised by the braid word, regardless of the nature of the fluid flow, i.e. if it is compressible, 
incompressible, viscous or inviscid. For three stirrers, the magnitude of the largest eigenvalue of 
the associated braid matrix, called the spectral radius gives a rate at which at least one material 



line is stretched. Finn et al. (2003) computed the stretching in time of a finite material line 
element and showed that the chaotic region of the fluid flow is of the same extent as the region 
of fluid the stirrers are moved through. 
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Figure 1. (A) Braid repre- 
sentation in space and time de- 
picting the exchange of stirrers 
for a pseudo-Anosov stirring pro- 
tocol. Three stirring operations 
are shown, corresponding to the 
exchange of stirrers <r 2 = (2,3), 
then a -i = (1,2) and again 
stirrers <7 2 = (2,3). The ro- 
tation stirrers trace out are de- 
picted in (B) for a finite-order 
stirring protocol and (C) for the 
pseudo-Anosov stirring protocol 
depicted in (A). 



In the remainder of this paper we present numerical simulations of three laser stirrers in an 
atomic BEC tracing out the finite-order and pseudo-Anosov stirring schemes A and B, presented 
above. Note the stirrers trace out a figure eight path, depicted in Fig. [2j We choose the velocity 
of the laser stirrers to be larger than the critical velocity for vortex nucleation in order to 
determine if a pseudo-Anosov stirring protocol (B) also results in a more random distribution 
of vortices than a finite-order stirring scheme (A). 
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Figure 2. Stirrer paths for 
first two stirring operations of 
A) finite-order stirring protocol 
and B) pseudo-Anosov stirring 
protocol. 



3. Modelling vortex dynamics 

The second quantized hamiltonian describing a trapped, weakly-interacting atomic Bose- 
Einstein condensate takes the form 



H 



df 



h 2 u 
-¥ (f ) — v 2 ^ (f) + (v (f ) + Vp) ¥ (?) * (f ) + — ¥ (?) ¥ (?) * (f ) * (f) 

2/n 2 



• (4) 



Here \fr (?) ' and \& (?) are the creation and annihilation operators for bosons of mass m, trapped 
in a potential V(P). The laser stirrers are described by Vp. The boson-boson interaction is 
given by a contact potential Uo = Airh 2 a/m, where a is the s-wave scattering length. The 
evolution of the field operator can be obtained from the Heisenberg equations of motion. At 
low temperatures, we can describe the condensed state by a classical mean field = and 
excitations (<f>), decomposing the boson operators as = '3/ + (0). The dynamics of the mean- 
field condensate wavefunction are then given by the nonlinear Schrodinger equation 



ih- 



dt 



2m 



V 2 ^ (f) + (V (?) + V P ) (f) + Uo\V (f) (f 



(5) 



We consider a condensate trapped such that the axial trapping frequency is much greater than 
the radial trapping frequency, uj z 3> u> r . In this case the condensate is quasi-two dimensional. 
For our simulations we model the dimensionless form of the 2D non-linear Schrodinger equation, 



dip 
l 8t 



1 



V 2 ^ + — tp + V P ip + K 2d \tp\ 2 ip . 



(6) 



Here the 2D interaction strength takes the form K2d = 2\/27raN/a z , where N is the total number 
of atoms in the condensate. The 2D condensate wavefunction is written as if) = a r ^/VN and we 
have imposed the normalization f dx|^| 2 = 1. The radial (axial) trapping frequency, co r (u) z ), 
determines the radial (axial) harmonic oscillator length a r = \fh/(muj r ) (a z = y/h/(muj z )). We 
choose the scaling for lengths and time to be r/a r = r and tu r = t, where f and t are dimensional 
with units of meters and seconds respectively. We solve Eq. §6§ pseudo-spectrally, applying an 
adaptive 4-5th order Runge-Kutta method in time. We select a stirring strength Vq ~ 2.6/x, 
where \i is the chemical potential of the BEC. We choose k = 10400, which taking 23 Na atoms 
corresponds to parameters of u z = 2ir x 50Hz, uj r = 2tt x 5Hz, N = 2.2 x 10 6 and a = 2.75nm. 
Simulations are run on a 512 2 grid of spatial extent —20 to 20. Snap-shots of time-evolution of 
the condensate density profile is depicted in Fig. [3} 



3.1. Decomposition of the Kinetic Energy contribution 

For a time-independent potential, the total energy of the system is conserved within the 
nonlinear Schrodinger equation. Writing the mean-field condensate wavefunction in terms of 
the condensate density n = |^| 2 , and phase 9, as tp = ^/nex^ [iO], the kinetic energy component 
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Figure 3. Vortex Dynamics: Condensate density profiles after simulation times given along 
the bottom of the figure. Magenta + and cyan — symbols denote vortices with positive and 
negative winding respectively. Row A.) finite-order stirring protocol at times corresponding to 
braid- words 02, 0102, 020102, (01 02) 2 and 02 (01 02) 2 '. Row B.) corresponds to pseudo-Anosov 
stirring protocols at identical times, corresponding to braid words 02, 0-102, 020-102, (0-102) 2 
and 02 (0_i02) 2 . 



of the total energy can be decomposed into contributions from the quantum pressure and kinetic 



energy density (Nore et al. (1997)) 



E KE = |dx QlWf) = Jdx Q|VV^| 2 + ^v| 2 ) . (7) 

Here we have defined the velocity field v = V0. The kinetic energy density contribution can be 
further decomposed into a compressible component for which the density weighted velocity field 
T = y/nv, satisfies V x T c = and an incompressible, divergence free component V • = 0, 
such that T = T' + T c . The total compressible energy is a useful measure of the total sound in 
the system and the total incompressible kinetic energy gives the total kinetic energy of vortices 

Efe = \ /dx|Y*>f . (8) 



One can subsequently define compressible and incompressible kinetic energy spectrums, 
however due to the small size of trapped atomic condensates, only a short range of length scales 
is accessible in these systems, in comparison to other 2D turbulent systems. Consequently, 
evaluating scaling exponents from such spectrums is of limited use and the derived scaling laws 
should be interpreted with care. We refrain from such an analysis here. 




Figure 4. Incompressible 
and compressible kinetic en- 
ergy distribution for finite- 
order stirring given by blue 
and green lines respectively, 
and pseudo-Anosov stirring, 
given by red and yellow lines 
respectively. 



3.2. Statistical measures of clustering 



Motivated by the seminal work of Ripley (1976, 1977), White et al. (2012) developed statistical 



measures of clustering that are appropriate for measuring if vortices are clustered or more 
randomly distributed in atomic Bose-Einstein condensates. These measures have the advantage 
of being experimentally accessible as they rely only on knowledge of the position and winding 
of vortices, taking into account both the sign and distance between a reference vortex i and its 
neighboring vortices j, up to its J th nearest-neighbor, 



Bj(t) 



1 bij (t) 

i=i j=i 



N 



J 



(9) 



N denotes all vortices that are at a distance di > Re from the condensate edge. If vortex i and 
its j th nearest neighbor are of opposite sign, bij = 0. If vortex i and its j th nearest neighbor 
are of the same sign and separation distance between them, dij < R c , then bij = 1. If vortex i 
and vortex j are separated by a distance larger than R c , (d^ > R c ), then b^ = 0. In this work 
we take the distance R c = Re to be greater than the average inter- vortex separation distance 
and of the order of the largest cluster size. These constraints are included in the statistical 
measure of vortex distribution in Eq. ^ , to account for edge effects that arise as we are dealing 
with a small system of finite size. The parameter limiting separation distance between vortices 
i and j gives a means to investigate the variation of clustering over spatial regions. We take 
vortex distribution to be clustered for values of Bj(t) > 0.5, while for more random vortex 
distributions Bj{t) < 0.5. Related measures of clustering have been recently applied to look at 
two dimensional turbulence in Bradley & Anderson (2012) and Reeves et al. (20126). 



4. Discussion and Conclusions 

We employ Eq. ( [9]) up to fourth-nearest neighbor vortices and track the corresponding statistical 
measure of clustering, which is plotted in Fig. [5} Evidently, neither the finite-order nor the 
pseudo-Anosov stirring protocol that we have chosen results in a vortex distribution that is 
significantly clustered. After a few stirring operations, the vortices appear randomly distributed, 
independent of stirring protocol. This result is an indication of the potential nature of quantum 
fluids. Although material line elements are expected to be exponentially stretched for pseudo- 
Anosov stirring in contrast to the linear stretching of material line elements in the finite-order 
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Figure 5. Evolution 
of vortex distribution statis- 
tics given by fourth nearest 
neighbor vortices for finite- 
order (blue *) and pseudo- 
Anosov (green *) stirring 
schemes. 
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case, this does not appear to influence the vortex distribution. In other words, the vortex 
distribution is independent of if topological chaos is built into the flow by the method of stirring. 
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Figure 6. Evolution of 
number of vortices with pos- 
itive and negative winding 
given by red diamonds and 
magenta circles respectively 
for finite-order stirring and 
blue diamonds and blue cir- 
cles respectively for pseudo- 
Anosov stirring. 
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From the analysis of the incompressible and compressible kinetic energy distribution (see 
Fig. [4]), it is found that pseudo-Anosov stirring results in less sound being generated in 
comparison to finite-order stirring protocols. The finite-order stirring scheme creates more sound 
in the condensate as result of the abrupt change in stirring direction during sequential stirring 
operations in this scheme. Although stirrers in both schemes follow trajectories that trace out a 
figure of eight shape, such abrupt changes in stirring direction do not occur in the pseudo-Anosov 
stirring scheme employed. Although at times t < 30 the total incompressible kinetic energy is 
initially similar for both stirring protocols, at later times more incompressible kinetic energy is 
generated for finite-order stirring. This difference is also reflected in the total vortex number for 
both stirring protocols (refer to Fig. [6]). Finite-order stirring generates a larger imbalance in the 
total number of positive and negative vortices than pseudo-Anosov stirring. The oscillations in 
the vortex number and kinetic energy distribution for finite order stirring reflect the oscillations 
in the total condensate volume induced by this stirring protocol. 



In conclusion, we have applied both pseudo-Anosov and finite-order stirring protocols to 
introduce vortices in a quasi-two dimensional Bose-Einstein condensate. Efficient mixing, as 
described by exponential stretching of material lines, could be thought to also translate into 
efficient vortex mixing. However, we observe that for a trapped inhomogeneous Bose-Einstein 
condensate, there is no appreciable increase in 'mixing' of vortices, as determined by the random, 
or clustered nature of the vortex distribution. That is, we do not find any significant correlation 
between if topological chaos is being introduced into the flow and how random the resulting 
vortex distribution is. In future work, we set up a stirring scheme that induces topological chaos 
in three dimensional quantum fluid flow and also analyze the onset of chaos for both finite-order 
and pseudo-Anosov stirring schemes. 
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